Abstract. Retrograde waves with frequencies much lower 
than the rotation frequency become trapped in the solar 
radiative interior. The cigcnfunctions of the compressible, 
nonadiabatic, Rossby-like modes (e-mechanism and radia- 
tive losses taken into account) are obtained by an asymp- 
totic method assuming a very small latitudinal gradient of 
rotation, without an arbitrary choice of other free parame- 
ters. An integral dispersion relation for the complex eigen- 
frequencies is derived as a solution of the boundary value 
problem. The discovered resonant cavity modes (called 
] .R-modes) are fundamentally different from the known r- 
f^i ■ modes: their frequencies are functions of the solar inte- 
' rior structure, and the reason for their existence is not re- 
(N lated to geometrical effects. The most unstable i?-modes 
£^ \ are those with periods of «l-3yr, 18-30yr, and 1500- 
^ ■ 20000yr; these three separate period ranges are known 
| from solar and geophysical data. The growing times of 

■ those modes which are unstable with respect to the e- 
! \ mechanism are ss 10 2 , 10 3 , and 10 5 years, respectively. The 

■ amplitudes of the i?-modes are growing towards the cen- 
f^*) | ter of the Sun. We discuss some prospects to develop the 
C^l ■ theory of i?-modes as a driver of the dynamics in the con- 

vective zone which could explain, e.g., observed short-term 
fluctuations of rotation, a control of the solar magnetic cy- 
' cle, and abrupt changes of terrestrial climate in the past. 
P : 

Key words: hydrodynamics - Sun: activity - Sun: inte- 
~ q . rior - Sun: oscillations - Sun: rotation 
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1. Introduction 

The 22-year magnetic cycle of solar activity is the most 
prominent phenomenon of several large-scale dynamic 
events that occur in the Sun. (Really, the magnetic half cy- 
cles or sunspot number cycles vary in length between 9-13 
years, and 11 yr is an average of the ~ 20 half-cycles avail- 
able.) An explanation of the basic mechanism underlying 
this famcous phenomenon is the fundamental challenge 
of solar physics. The achievements of the theory of the 
a-u dynamo turned out to be a great success. However, 
neither all observations of magnetic and flow fields nor 
the radiation fluxes which are related to this phenomenon 
and which are measured at the surface of the Sun or indi- 
rectly, by helioseismology, in its interior, can be explained 
unambigously in this way. Although our present work is 
not directly related to the dynamo theory, we will outline 
here those difficulties which have common points with our 
results. 

1.1. Some problems of dynamo theory 

As a consequence of our imperfect knowledge of basic 
characteristics of turbulent convection as well as merid- 
ional circulation and details of the rotation of the Sun's in- 
terior, the solutions of the dynamo equations become func- 
tions of many free, unknown parameters (e.g. Stix 1976). 
For instance, by clever combinations of these parameters 
it is possible to get from kinematic theory an oscillatory 
magnetic field with a 22-year period and a growing ampli- 
tude. However, another choice of these parameters leads 
to waves of growing amplitude for other periods. So one 
could draw a butterfly diagram not only with an 11-year 
periodicity. It remains still an open question which of the 
clever combinations resulting in a solar-like 22-year activ- 
ity cycle is realized in the Sun. We could not find a work 
on the dynamo wave problem, showing that just the 22- 
year period is preferred among others with a maximum 
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growth rate and with the spatial scales required for solar 
activity. Instead, many authors pointed out that the cycle 
period of 22 years is hard to explain (Stix 1991; Oilman 
1992; Levy 1992; Schmitt 1993; Brandenburg 1994; Weiss 
1994; Riidiger & Arlt 2000). 

From the solution of the inverse problem of helioseis- 
mology (e.g. Tomczyk et al. 1995) it is known that the 
convective envelope of the Sun is rotating with a latitude 
dependence of the angular velocity similar to that of the 
surface but almost rigid in radial direction. A stronger ra- 
dial gradient which is required for the a-uo dynamo mech- 
anism is located in a shallow layer (thickness « 0.05i? Q 
(Kosovichev 1996), where R@ is the solar radius) immedi- 
ately below the convective zone — the tachocline (Spiegel 
& Zahn 1992). Below the tachocline up to a depth of 
at least O.5i? the radiative interior is rotating with an 
angular velocity law similar to that of a solid-body. The 
question arises: what compels the Sun to rotate in such a 
strange manner, which is different from the generally ac- 
cepted, theoretically predicted stable rotation law? How 
to handle a dynamo theory for which the 'u>' area is sep- 
arated from the l a' area over a large part of the extent of 
the convective zone? In order to solve this problem Parker 
(1993) has put forward the idea of an interface dynamo, 
the basic features of which existed already in earlier dy- 
namo models (Steenbeck et al. 1966). To close the cycle 
of such a stretched dynamo it is necessary to have some 
mechanism delivering toroidal magnetic flux, arising by 
the shear of differential axisymmetric rotation (Cowling 
1953) in the tachocline, to the 'a' dynamo area (e.g. Mof- 
fatt 1978; Krause & Radler 1980). To get a solar-like mag- 
netic activity it is necessary to suppose the existence of a 
huge 10 5 G) toroidal magnetic field to create enough 
magnetic buouancy for the leakage of magnetic flux and 
to solve the tilt problem of lifting loops (e.g. Caligari et 
al. 1998). Moreover, a high magnetic diffusivity contrast 
between the convective envelope and the underlying ra- 
diative core should be assumed to solve the quenching 
problem of the a effect (see, e.g., Fan et al. 1993; Catta- 



N.S. Dzhalilov, J. Staude & V.N. Oraevsky: Long-period eigenoscillations of the Sun. I. 



3 



neo & Hughes 1996). However, it is a major challenge for 
any dynamo model to produce such strong fields. 

The idea of the interface dynamo was further devel- 
oped, e.g. by Charbonneau & MacGregor (1997). Later, 
a fit to the real solar rotation profile with its latitudinal 
and radial dependencies has been included by Markiel & 
Thomas (1999), but so far no satisfactory solar-like oscil- 
latory solutions for the interface dynamo have been found. 
Growing wave solutions are suppressed by the latitudinal 
shear. 

1.2. Spinning- down rotation problem 

Mechanisms for braking the solar internal rotation are 
also under discussion. The character of the core rotation is 
not clear because here the accuracy of hclioseismic inver- 
sions gets worse (Chaplin et al. 1999) and the results seem 
to be in contradiction with the oblateness measurements 
(Paterno et al. 1996). There are some suggestions that a 
deceleration of the radiative interior depends on the trans- 
port of angular momentum between this region and the 
convective zone. For instance, Mestel & Weiss (1987) sup- 
posed that even a weak large-scale magnetic field would 
be sufficient to couple very efficiently the interior and the 
convective zone, leading essentially to solid body rotation. 
In this way the magnetic torques can also extract angular 
momentum from the radiative interior (e.g. Charbonneau 
& MacGregor 1993). 

The wave mechanism for the solution of this problem is 
more popular. Schatzman (1993), Zahn et al. (1997), and 
Kumar & Quataert (1997) have concluded that the solid 
rotation of the radiative interior is a direct consequence 
of the effect of internal gravity waves. Gravity waves gen- 
erated near the interface between the convective and ra- 
diative regions transport retrograde angular momentum 
into the interior, thereby spinning it down. Here the main 
idea is that the isotropically generated gravity waves be- 
come anisotropic due to Dopplcr shifts of frequencies in 
the differentially rotating Sun. In that way for anisotropic 
retrograde and prograde waves the radiative damping is 
different, and the residual negative angular wave momen- 
tum may compel the solar radiative interior to co-rotate 
with the convective zone. This idea has been further de- 
veloped by Kumar et al. (1999) including a toroidal mag- 
netic field to explain the existence of the unstable shear 
layer 'tachocline'. However, Ringot (1998) has shown that 
a quasi-solid rotation of the radiative interior cannot be a 
direct consequence of the action of internal gravity waves 
produced in the convective zone. Cough (1997) questioned 
this idea emphasizing that the mechanism can work only if 
the waves are generated with strong amplitudes to trans- 
port the required angular momentum. This means, reso- 
nance waves are required, but such waves may penetrate 
only to distances less than lO _5 i?0 beneath the convective 
zone due to the strong radiative damping. These waves 



must deposit their angular momentum before returning 
to the convective zone, but not before penetrating far into 
the radiative interior. 

For the wave mechanism the question of an anisotropic 
propagation relative to the azimuthal rotation is a key 
moment. Fritts et al. (1998) have shown that convec- 
tion, penetrating into the stratified and strongly sheared 
tachocline, can produce preferentially propagating gravity 
waves. 

There have also been speculations that the rotation of 
the core may be variable, perhaps with a time scale of the 
solar cycle (e.g. Cough 1985). The present paper is along 
these lines. 

From our short discussion we conclude that the con- 
vective envelope and the radiative interior are coupled to 
each other through a certain global agent, resulting in al- 
most co-rotation. To advance the solution of the problem 
the dynamo theory should take into account the presence 
of this global agent. We suppose that really this agent is 
provided by waves with the following properties: 

Waves should represent large-scale global eigenoscilla- 
tions of the Sun. Their origin must be related to rotation, 
they must be strongly anisotropic with respect to the az- 
imuthal angle. Looking at the characteristics of the so- 
lar cycle we immediately see the high coherency of these 
global motions (the constant periods, phase shifts, ampli- 
tudes, the latitude appearence, etc.). Activity grows in the 
first phase with a timescale which is considerably shorter 
than the decay time in the second phase; this fact and 
the quick eruptive release of energy by the reconnection 
mechanism indicate that the waves must be unstable. 

It is noteworthy that the inner gravity waves do not 
fulfill these requirements. The quesion is whether the r- 
modes do? 

1.3. r-modes 

In a non-rotating sphere (il = 0, where Q is the angu- 
lar frequency of solar rotation) the wave motion is subdi- 
vided into two non-coupling components: spheroidal p— , 
/— and g— modes (for which the main restoring forces 
are pressure gradient and buouancy) and toroidal modes 
(e.g. Unno et al. 1989). Toroidal modes are degener- 
ated horizontal eddy motions confined to a spherical sur- 
face with a radius r for which u = 0, divv = 0, and 
v = Qp(r) x (0,^- 9 ^,-^)Yr(e,4>). Here Y™ is the 
spherical harmonic with a degree I and order m, 9 is the 
colatitude, <f> is the azimuthal angle in the spherical polar 
coordinates, v is the fluid velocity field, w is the angular 
frequency of the fluid motion, and Qf 1 ^) is an arbitrary 
amplitude function. Toroidal modes have zero radial ve- 
locity but have non-zero radial vorticity, (rotv) r ^ (for 
the spheroidal modes it is vice versa). These modes do not 
alter the equilibrium configuration. 
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When a slow rotation (fi 2 < f2 2 = GM/R?) is included 
the spheroidal modes are slightly modified but they keep 
their main properties. Degeneracy of toroidal modes is re- 
moved only partially by the rotation, and quasi-toroidal 
waves - known as r-modes - appear with a non-zero fre- 
quency of a; ss 2Qm/l(l + 1) in the rotating frame (Pa- 
paloizou & Pringle 1978; Brayn 1889). Usually the govern- 
ing equations of the r-modes are obtained by expanding 
the initial physical variables of the equations in the ro- 
tating system into power series with respect to the small 
parameter (f2/f2 g ) 2 (s3 1CT 4 ' 7 for the Sun, e.g. Papaloizou 
& Pringle 1978; Provost et al. 1981; Smeyers et al. 1981; 
Saio 1982). These power series practically describe the de- 
viation of the surface of the star from its initial spher- 
ical state, resulting from rotation through Coriolis and 
centrifugal forces. As a result of the deformation of the 
spherical surface with a radius r the radial vorticity of 
the toroidal modes cause a surface pressure perturbation 
through the Coriolis force. However, the r-modes practi- 
cally keep the main properties of toroidal flows: « r k 0, 
divv ss 0. The degeneracy of the r-modes is that their 
frequencies hardly depend on Q" l (r), i.e. they are inde- 
pendent from the inner structure of the star. For the I = 1 
modes the frequency in the inertial system is again close 
to zero, lo ss (Papaloizou & Pringle 1978). The r-mode 
equations define the amplitudes Q[™(r), and taking into 
account the next terms with small VL 2 /SI 2 in the series 
practically does not change the frequencies. 

Due to the fact that the r-modes are practically surface 
deformation waves, some similarity of these waves to the 
surface gravity waves or to the /-modes is apparent. For 
high I the /-modes are an analogy to surface gravity waves 
in a plane-parallel fluid with lo 2 = gk±. In the Cowling ap- 
proximation /-modes with I = 1 have zero frequency too, 
lo ss (Unno et al. 1989). This corresponds to a parallel 
displacement of the whole star. For high I the /-mode fre- 
quencies are also independent of the inner structure, with 
lo 2 ss IVL 2 (Gough 1980). So, r-modes are also fundamental 
rotating modes with an inertial frequency, lo < 2S7. 

For the Sun the properties of r-modes have been in- 
vestigated in great detail by Wolff et al. (1986) and Wolff 
(1998; 2000; and refs. therein). 

Some properties of the r-modes are also similar to 
those of the Rossby waves in geophysics (Pedlosky 1982). 
Similar to the Rossby waves and unlike the g-modes the 
r-modes are strongly anisotropic. They propagate only in 
azimuthal direction, opposite to rotation (i.e. they are ret- 
rograde waves in the co-rotating frame). Because we are 
interested in length scales corresponding to those of large 
sunspots, we have to consider r-modes with / ss 100. To 
get oscillations with periods of years (lo/Q ss 10~ 2 ) wc 
must choose m ss I 1, just such r-modes are physically 
more interesting (Lockitch & Friedman 1999). However, 
in the case of high I the amplitudes of the r-modes will 
be concentrated near the surface of the Sun (Provost et 
al. 1981; Wolff 2000), and so they can actively interact 



with convective motions (Wolff 1997; 2000). Because for 
these modes v r ss and div v ss 0, their chance to take 
part in the redistribution of angular rotation momentum 
in the radiative interior is low. Note that the slow solar 
differential rotation does not change the behavior of such 
r-modes with m = I > 1 (Wolff 1998). 

Looking for further analogies between waves connected 
with gravity and with rotation, we remember that beside 
the surface gravity waves there exist internal gravity waves 
with lo 2 ss N 2 k 2 /(k 2 + k 2 ), the frequencies of which de- 
pend on the inner structure (N is the Brunt- Vaisala fre- 
quency). Similar to these waves there exist 'true' Rossby 
(not deformation) waves, the frequency of which depends 
also on the internal structure. 

1.4- Rossby waves 

We include here a short review on the main features of 
Rossby waves; they have been investigated in great de- 
tail in geophysics (e.g. Pedlosky 1982; Gill 1982). In the 
simplest case, that is in a plane-parallel, homogeneous, ro- 
tating layer, the dispersion relation for the Rossby waves 
is lo ss 2Vtj3k x /{k 2 + k 2 + k 2 ). Here k x is the wave number 
perpendicular to the rotation axis, k z is expressed by the 
internal deformation radius of Rossby which depends on 
the Brunt- Vaisala frequency, /3 is the transverse gradient 
(in y direction) of the Coriolis parameter: a vertical com- 
ponent of the 'planetary' vorticity 21~2 in the given local 
point. Unlike the r-modes the Rossby wave frequencies 
are functions of the internal structure and have maximum 
dependence on the gradient (3: lo — > if k x — > and if 
k x — > oo. 

Any disturbance of the local flow in a rotating frame 
may generate waves of the Rossby type. These waves exist 
only if there is a gradient of the potential vorticity II = 
(u) a ■ VH)/p. Here an absolute vorticity is the sum of the 
relative and the planetary vorticities, to a = rotv + 217, 5 
is any conserved scalar quantity, ^ = (for instance, for 
adiabatic motion that could be the entropy or the density 
in the case of incompressible plasma). The Rossby wave 
motion is a solution of the nonlinear equation for transport 
of n. The potential vorticity is conserved if the medium 
is barotropic (Vp x Vp = 0) and if there are no torques. 
The rotation of the frame is added to any vorticity in the 
velocity field. Any motion within a rotating fluid serves as 
a potential source for vorticity. 

The relative vorticity may be evoked by the geomet- 
rical surface as well as by internal gradients. It depends 
on the choice of the function S(r, 6, <j>) and on J7(r, 0). For 
example, an unevenness of the ocean bottom causes the 
topographic Rossby waves, or a dependence of the Coriolis 
parameter on the earth latitude (F — 2f2 sin ip, where ip is 
the geographic latitude) is the main cause of atmospheric 
Rossby waves. 

In the solar dynamo context the ability of Rossby 
waves to induce solar-like magnetic fields has been consid- 
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ered by Gilman (1969). Here the mechanism for sustain- 
ing the Rossby waves is a latitudinal temperature gradi- 
ent in a thin, rotating, incompressible convective zone. To 
interpret the dynamical features of large-scale magnetic 
fields the Rossby vorticies excited within a thin layer be- 
neath the convective zone are considered by Tikhomolov 
& Mordvinov (1996) as the result of a deformation of the 
lower boundary of the convective zone. 



1.5. R-modes 

From the discussion in Subsection 1.4 we conclude that 
just Rossby-like waves could be suitable for our require- 
ments. As the main driving mechanism we choose a latitu- 
dinal (or horizontal) differential rotation, O = 0(6). Baker 
& Kippenhahn (1959) have pointed out that the uniform 
rotation of a star is not a typical case. Low frequencies 
(periods of years) could easily be obtained searching for 
the eigenoscillations of the Sun's radiative interior, where 
the gradient of the rotation speed is close to zero (in accor- 
dance with the hclioseismology results) . Large scales such 
as those associated with sunspots (k x R Q ~ 100) decrease 
the frequencies too. Similar to the r-modes the Rossby 
waves are strongly anisotropic (retrograde waves), but un- 
like the r-modes these waves are concentrated close to the 
solar center. These results have already been obtained by 
Oraevsky & Dzhalilov (1997), who investigated the trap- 
ping of adiabatic, incompressible Rossby-like waves in the 
solar interior. In the present work we take into account 
compressibility for the nonadiabatic waves. We look for 
unstable waves. It is clear that the necessary condition for 
the Kelvin-Hclmholtz shear instability, AN 2 — (r^) 2 < 
(Ando 1985), is not fulfilled in the radiative interior. Then 
we decided to include the thermal e-mechanism of insta- 
bility which is favoured at low frequencies (Unno et al. 
1989). To balance the e-mechanism the radiative losses in 
the diffusion regime are included. To exclude all geometri- 
cal effects we ignore the influence of the spherical surface 
at the given radius. The dispersion relation in the limit of 
adiabatic incompressiblity and at very low frequencies is 
the same as that for Rossby waves in geophysics. In or- 
der to distinguish these rotational body waves from the 
r-modes we call them R-modes (Rossby rotation). 

The governing fourth order equation is obtained from 
the basic equations in Sect. 2. Some qualitative analysis of 
the wave cavity trapping is done for the simpler adiabatic 
case in Sect. 3. Using the asymptotic solutions obtained 
in Sect. 4 the complex boundary value problem is solved 
in Sect. 5. The calculation of the eigenfrequencies and the 
instability analysis are done in Sect. 6. The obtained un- 
stable modes are shortly discussed in Sect. 7. 



2. Setting the problem 



2.1. Basic equations 

Let us investigate global motions with large timescales 
such that the Rossby number is small, w/fi < 1. Before 
the appearence of any disturbances the basic stationary 
state of the rotating star is defined mainly by the bal- 
ance of pressure gradient, gravity force, and forces ex- 
erted by the noninertiality of the motion (Coriolis and 
centrifugal forces). In the case of an incompressible fluid 
with a homogeneous rotation rate usually this state is 
called 'geostrophic balance'. A star disturbed by an ex- 
ternal force tends to return to this basic state. Our aim 
is to study for the Sun the dynamics of small deviations 
from the steady geostrophic balance. For this purpose it 
is natural to write the dynamic equations in a frame ro- 
tating together with the Sun. The magnetic field will be 
ignored. For arbitrary 0(r, $) the equation of momentum 
in conventional definitions is given by 

dv „ 1„ dfl „ , „ , . 

- — \-2flxv = — Vp+g+rx — — fix (f2xr)+/i„V 2 'w,(l) 
dt p dt 

where p, v is the kinematic/turbulent viscosity coefficient. 
In the next steps this equation will be simplified keep- 
ing the main features of the motion. We consider linear 
waves without taking into account convective and merid- 
ional flows, Vo = 0, and we suppose that the angular ve- 
locity O does not depend on time. For the basic state we 
have from Eq. ([!]) 

— — Vpo + g = n x (n x r). 

Po 

However, we can exclude practically everywhere in the Sun 
centrifugal acceleration and consider it as a small correc- 
tion to g. Then the spherically symmetric basic state is 
defined by Vpo ~ Po<7- For the linear oscillations we have 
from Eq. (Q): 
dv 

— + 20 x v + r sin • Vfi)e v = 



dt 



Po 



-Vp' 



Vp , 
Po 



p, v V 2 v. 



(2) 



Here e v is the unit vector in azimuthal direction, and 
variables with a prime are Eulerian perturbations. Eq. (J2J) 
is written in rotating spherical polar coordinates (r, <&, tp). 
It coincides with the equation of motion by Unno et al. 
(1989) which is derived for an inertial frame, if the oper- 
ator d/dt is replaced by d/dt + fid /dip. To simplify our 
further discussion we shall use the Cowling approxima- 
tion, g' = 0, which has a sufficient degree of accuracy for 
the analysis of short waves. 

The next simplification of Eq. (Q) is connected with 
the quasi-rigid rotation of the inner part of the Sun below 
the convective zone, which is known from the solution of 
the inverse problem of helioseismology. In this case we 
can omit the third term of the l.h.s of Eq. (||). Such a 
restriction of the gradients of fi(r, i?) requires to obey the 
conditions: 

dO, dO, 
2fi>r— , 2ft»tan(tf) — . (3) 



dr 
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To obtain these conditions we used the components of the 
vector ft in spherical polar coordinates, J~i = {i7cosi9, 
-fisini?,0}. 

The next equations are the mass and energy conserva- 
tion equations in the standard form: 



dp 

— + pdwv = 0, 
at 

d T 

CvP—rr+pCLiYV = -£, 
at 



(4) 
(5) 



-(7-1)/;. 



where: - C = Q{p, T) - V • q R - V • (K T 5/2 VT). 
Another form of Eq. (||) is 
dp _ 2 dp 
dt C " dt 

Here 7 = c p /c v , c p and c v are the specific heats at constant 
pressure and volume, respectively, c s is the sound speed, 
the source function Q(p,T) is the sum of nuclear and vis- 
cous heat generation rates per unit volume: p(e^ + £ v ), 
qR is the radiative energy flux, KqT 5 / 2 is the coefficient 
of electron heat conductivity. In the interior of the Sun we 
shall neglect viscous heating (e v — 0). For the power of 
nuclear reactions we have Q ~ p 2 T a . In particular for the 
p-p reaction Q = pe pp w 9 x 10 _30 x#p 2 T 4 , where e pp is 
given in [erg/gs] and xh = 0.73. 

In the limit of incompressible fluid, dp/dt — or 
divu = 0, it follows from Eq. (§) that the condition 
c 2 — > 00 is not needed to satisfy dp/dt ^ 0. Hence, in 
a dissipative (C 7^ 0), incompressible fluid sound can- 
not propagate instantaneously. It means, we cannot use the 



condition of 



00 to get the incompressible limit for 



nonadiabatic waves . 

Now we will try to simplify the energy loss function 
C assuming reasonable approximations for the Sun's in- 
terior. We use the formula for the heat conductivity of a 
fully ionized gas to show that in the Sun radiative trans- 
port of energy is more important than that by particle 
heat conductivity: 

k t = K Q T 5/2 w Kr 6 T 5/2 erg/scmK. 

The radiative flux is given by the radiative diffusion equa- 
tion 



Qr = -K fi VT, kr 



3\ 



(J) 



where kr is the radiative heat conductivity, as = orc/A = 
5.67 x 1CP 5 erg/cm 2 K 4 s is the Stefan-Boltzmann con- 
stant, and x is the Rosseland mean opacity. For the 
opacity we will use Kramers' formula: x = 3.68 x 
lQ 22 p 2 T~ 7 / 2 cm -1 . Then, for the radiative conductivity 
we have 



.22 x 10 



j.13/2 

' 7 — k— erg/scmK. 



Now we find that k r /k t ~ 10" 20 T 4 /p 2 > 10 3 . This con- 
dition is fulfilled in the whole Sun. That means, we can 
exclude the heat conductivity term from Eq. (p|) . The next 
simplification is connected with the ratio T*l 2 j p which 
is practically not changed over the radius and « 10 16 . 



It means, we can introduce the constant ae = 8.22 x 
\Q- 27 T b /p 2 fa 8.2 x 10 5 . Then we have 
kr Ri asoT 3 / 2 erg/scmK. 

Now for the right-hand side of Eq. (|^) we have 
-£^ae V 2 T 5 / 2 + Q. (8) 

The main non-perturbed energetic state of the Sun is de- 
fined by the condition Cq = 0. In our case this condition 
is 



d 2 T, 







Qo = ~^o dr2 



(9) 



(6) 2.2. Choice of the frame of reference 



In the next step we try to get the analytical solutions of the 
wave equations and to solve the boundary value problem. 
For this aim we shall investigate the short waves (WKB) 
approximation for which the effects of curvature are unim- 
portant. It means we can apply the rotating plane-parallel 
stratification approximation. Such an approach has its ad- 
vantages and disadvantages. The main advantage is that 
at the end points of the integration path (center and pole 
of the Sun) we have no singularity, and this gives a chance 
to find the solutions analytically. The disadvantages are 
connected with the following: a) long waves are excluded 
but they are very important, for instance, for the transfer 
of angular momentum; b) the lost 'end' singularities cor- 
respond to a real physical behavior of the waves; c) in this 
approach we get two distinct directions, £1 and g (taking 
into account all components of ft), which are absent in a 
self-gravitating, radially stratified sphere. The deviation 
of the direction of stratification of the plane-parallel fluid 
layer from the direction of its axis of rotation should lead 
to additional, physically doubtful results. In geophysical 
hydrodynamics this problem is solved by applying the i /3- 
plane' within the frame of the traditional approximation 
(Pedlosky 1982; Shore 1992), where the component of O 
parallel to g in the given surface point (local vertical) is 
retained only in the governing equations. Here we will use 
the same approach. 

Let us take an arbitrary point at the surface of the 
rotating sphere. The position of this point is defined by 
its radius r, its co-latitude 1?, and its azimuth angle. We 
assign to this point a local, left-handed Cartesian system 
of coordinates {x, y, z}, where the z-axis is directed along 
the radius (local vertical), the direction of the y-axis is 
meridional (towards the pole), and that of the cc-axis az- 
imuthal. In this frame of reference Jl = {OjO^fiz} = 
{0, ft sin$, f2 cos Strictly speaking, the z-axis coincides 
with the rotation axis only at the pole ($ = 0). In the case 
of a homogeneous fluid Q is included into the wave equa- 
tion only in the term 

d 2 d 2 ^ d 2 



(nvy = n 2 — + 2n y n z 



Here we can neglect Vt y (O y = 0), if the condition 
<C Mz;|j| is fulfilled. This condition is named 'tra- 
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ditional approximation'. For harmonic motions we have 
-jjy- = ik y , = ik z , and the traditional approximation 
corresponds to \k y \ <C \k z \ cottf. If the spatial scale of the 
wave motion in vertical direction is much smaller than in 
horizontal direction (at latitudes not close to the equator) , 
we can restrict ourselves to retain only Q z in the governing 
equations. This condition for the traditional approxima- 
tion remains valid, if a radial stratification is included (Lee 
& Saio 1997). 

In order to construct the '/3-plane' limit we expand 
Q z around a fixed ?9 = #o + 8ti): Q z = 0(??)cosi? w 
(1 + (3y)£l cosiV Here y = R e Stf, fl = O(tf ), and 

- tani9 ) — . 

Hq 



n &d 

Rossby waves are possible only if p ^ 0. The parameter 
P is a sum of two terms: the second one (tani?) is due 
to the geometrical change of the Coriolis parameter with 
latitude. This term exists always, even if the rotation is 
rigid. The r-modes are connected with this term. The first 
term in P appears if the differential rotation is considered, 
and the i?-modes are connected with this term. Close to 
the pole (t?o ~* 0) the first term is dominant. 

Note that such a '/3'-limit is applicable also around 
the equator plane, where the traditional approximation 
does not fit. The advantage of this limit is the possibil- 
ity to include $ as a parameter in the Cartesian system. 
In this way we use here an inertial Cartesian system of 
coordinates (x,y,z) in a frame rotating with an angular 
frequency Q(y,z). All non-perturbed model variables are 
functions of z only, and for the gravity acceleration we 
have g = {0, 0, — g(z)}. For the observer from the non- 
rotating frame the elements of fluid are moving due to ro- 
tation with a velocity Vo = ft x r = {— fly, fix, 0}, where 
in the frame of our approximation f2 ss Q z . Vq x < means 
that the x-axis is directed opposite to the rotation. 

2.3. Oscillation equations 

For linearization each physical variable / = /o + /' is 
decomposed into a mean term fo and a small fluctuating 
term /'. Neglecting terms of higher order than the first 
one we get our oscillation equations 



dp' ^ dp 
dt z dz 



po divv = 0, 



dv 
~dt 



2Slxv 



Po Po 



p v V v, 



dn 



dt 



dz 



CvPO I r ^T + V z ^- ) + P0 divt) = 33 V 2 ( Tq /2 ^— 



V 



Tn 



Qo[2^ + a^) =-£'. 
Po To, 1 



(10) 
(11) 

(12) 



For adiabatic oscillations £ = 0. We approach this regime 
by setting formally 33o = (Eq. 9). But in the nonadi- 
abatic case, which is considered here, £' is the sum of 
radiative damping and the e-mechanism terms. Eqs. (|lfj|)- 



(13) 



( |l2"| ) must be completed then by the equation of state, 
p = p(p, T). For an ideal gas we have 

£ = pL + TL. 

Po Po To ' 
For the ideal gas po = R g poTo/ p, m , where R g is the gas 
constant, [i m is the molecular weight, and the squared 
adiabatic sound speed c? s = 7c 2 is defined by the squared 
isothermal sound speed c 2 = po/ Po- 
We define 

33 



Qo = 



83t 



Qo 
CvPoTo' 
J_dJh 
To dz ' 
1 dp _ 
Po dz 



33 



CvPoTo ' 
1 dp 
po dz ' 



1 

—33 

7 



p 



3C, 



33t 
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Here N 2 (z) is the squared Brunt- Vaisala, frequency. 

and 33t are the reciprocal pressure, density, and 
temperature inhomogeneity scales, respectively, and Qo is 
the reciprocal of the characteristic Kelvin-Hclmholtz time- 
scale — a deviation from the thermal balance of the star is 
restored during this time. As the coefficients of Eqs. (|Io|)- 
( |l3|) are independent of the time t and of the space coor- 
dinate x we can set 
d d 

_ = _ = lkx . (14) 

Now we exclude p' from the system of Eqs. (|!Q)-(ll3|) and 
have 



iuj 



P_ 

Po 



To 



33 fl V z 



0, 



Po 

-c 2 . 
dy po 

2 d p' 



20', 



v 



Dv x = 0, 



d 



zpo 



2il v x 

r 

Jo 



DVy 



■Dv, 



0, 



0, 



(iu> + 2Qo) \-ee p v z +^u = f, 



Po 



u = divv = ik x v x 



OVy 

dy 



dv z 
dz 



(15) 
(16) 
(17) 
(18) 
(19) 
(20) 



where the operators D and / are defined as 



D = iuj- 



f = (a- 



2)Qo ¥o 



aeoV 2 ^T 5/2 -j,(21) 

and the relation po/{c v poTo) = 7 — 1 is used. In the next 

steps the viscosity appears in the coefficients of the equa- 

2 

tions in the form + p, v ■ 0(1). For a fully ionized plasma 
the kinematic viscosity is p v w 10 _16 T 5 / 2 /p [cm 2 /s]. For 
the solar situation \x v w O(l), if we do not take into ac- 
count the turbulent viscosity. As we are interested in very 

2 

low frequencies the condition /i v is valid. Then we 

can put D ~ iu> (non- viscous case). Now Eqs. (16)-([l7|) 
define the horizontal components of velocity and hence its 
two-dimensional divergence div^v: 



T 



l -a 2 



1 



1 d 

ak x dy 



EL 
Po' 



(22) 
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oj I -a 2 
a 1 ~ p' 

divj_t; = c 2 ik x — M — . 

Pa 

„2 



where 

M-1 + - — 

k x dy 



k x dy) po ' 



dy 2 



(23) 

(24) 

da 
dy' 



k x <r> \dy 2 V 2n 
Excluding u = div±v + dv z /dz by using Eq. ( p"5[ ) we have 
the system of equations 

r2 d_i_ V 

dzpo 
d 

+ 7T 

oz 



iui v z 



9 



7fr - s(y,z)— + — 
T po iui 



To' 
-0, 



(25) 
(26) 
(27) 



2Qo 



(28) 



(29) 



(30) 



0, it follows immedi- q 



p' IN 2 T' 1 - 

Po iuj g T 7* 

where the dimensionless quantites are 

5=1- c 2 k x ^—M, 7* = 7«w, q = - 

uj 2 7 7* 

Now we use the condition a <^ 1 and receive 
c 2 d p' 

2£1 dy po ' 
_ i II 

Wy ~ 2n x p 

In the case of rigid rotation, § 

ately from Eqs. (P9|)-(|30|) that divj_t> = 0. It means, that 
the incompressible case (divu = 0) for which everywhere 
v z = (from the boundary condition at the center) is 
not so interesting for astrophysical situations. In this case 
p' = and the equation of heat conductivity, Eq. ((|l2|), 
is separated from the equation of motion. In geophysical 
situations just this case is interesting, when geostrophic 
eddies are investigated. For our task waves with v z =/= 
are more important. 

Our next step is to separate the z and y dependence 
of the variables in the governing equations to have finally 
one ordinary differential equation. In the general case such 
a separation is not possible, and we consider here the very 
simple case when the function g^(jj) is independent of 
y. Only in this case it is possible to separate the equa- 
tions with respect to the variables y and z and we can 
write Jp- = — k 2 . Here k y should be determined from the 
boundary conditions, and a complex k y is not excluded. In 
this way the system of partial differential equations in the 
plane-parallel approximation is reduced to ordinary dif- 
ferential equations. Now we assume the following formula 
for the rotation profile 



(31) 



1 + (3y/R Q ' 

where (3 > 0. This profile reveals that the rotation rate 
decreases at the pole. If y — > oo then fi — * 0. Helioseismol- 
ogy inversions predict f2(i?) in the convective zone which is 
almost the same as that at the surface (/3 « 0(1)), but be- 
low the bottom of the convective zone (at the tachocline) 
rotation is close to a solid-body rotation (/? <C 1). As 
y» = y/i?© < 1, in the deepest layers of the Sun we can 



assume the rotation rate Q(y,z) « J7(z)(l — /3y*), where 
n(z) sa O(n ). Then we have in Eq. (28) 

M « M = 1 + , ^| x , S = s » 1 - M^|4, (32) 

where kj_ = k 2 + k 2 . 

Our next step is to derive one differential equation for 
the temperature perturbations. The variable v z is easily 
excluded from Eqs. (|25|)-(|27|). Then we get for the pressure 
perturbations 

N 2 ' 



V 



Pa 



H 

To~ 



w q 



d_ 

dz 



To~ 



TL 
J_ 

7* 



7* 



(33) 



1 



d{z) 



where 



d_ 

dz 



N' 



^ 3. d = — = ^-»p-(34) 
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The parameter d is positive in the solar radiative interior, 
d > 0, while we have d < in the convective zone. The 
pressure perturbations (Eq. 33) have a singularity at a = 
0. However, it will be shown below that this singularity is 
removable. 

Introducing a new dependent variable 



K T~o' K 



T, 



on 



5/2 



To, 



const, 



we receive the final equation of fourth order 

e"" + a 3 e'" + a 2 e" + a x e' + a q e = o, 

where O' = d0/dz,e o = &oT^ 2 H* and 



Ao = aia 2 
A\ = aia 3 + a 2 + a' 3 , 
A 2 - 



£oC" 



AT2 



(35) 
(36) 

(36') 



a3 



ai6 2 - 6oj ^3 



\n(poT ad 2 



a' d' 

-ax = h — 

a a 

1 

a 2 = — 



h - -3ST 



qui 
dc 2 ' 



a-3 = 



— - k 2 

£o 



7 

"(2 



2 /7V2 



33l, 



r "^ 1 — 



9 



- 1 
1 



d! qoJ 2 

4 kT_ 

5 K 



5 k 7 

In order to solve Eq. (|3^) we normalize it to get a 
dimensionless equation. For this purpose we normalize the 
radial distance to the solar radius: z z/Rq and get 
finally the following equation 
e0"" + e^ 3 0'" + (l + e^2)0"+ 

+ (?/>i+£^i)0' + (V>o+ f^o)0 = O, (37) 
where we kept the notation k\ Rq^, the inverse scale 
heights ae p ,T are defined as above (below Eq. (13)) but 
with a normalized z, and 
^o = -aQ 2 (N 2 /cu 2 - 1 + ek 2 ) + bx{a x + ae T + b^/b^, 
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<Pi 

"¥>2 



dz 



¥>3 = ™^J ln (poTaad 2 ) 



ai(aei + b 2 k 2 ) + (ae'j + b 2 k 2 ± )', 
2se^, 



(37') 



4 k" 
-k 2 ± ip 3 + -ae p — 



- -1-2 



= k 2 , + a\b 2 + 6 2 — u) 2 a, 



-ax = — 



d' 

— + ae. 
a 



p 



~i1 

" d' q = 



7-1 



4 tf_ 

5 K 



aeT 



38t 



7 



= 1-A 



~2<7 

d' 



-9 
"d 



-1 



71/ = 1 



TV 2 
2~ 



1 



d_ 

dz 



£0 

k x ,±Re, 



ss T 



5/2 



1 



fc^Af 
ujRq/c, A — p^^- 



2ftwi? Q 

Remind that 7* = jilu. Here the main parameter is A, it 
includes the rotation rate gradient (/3=const). 

Because we are interested in very low frequency oscil- 
lations with periods of 1-20 years, we take u> s» 10 -8 s -1 , 



n 5/2 

Pa 



const = 10 



2 x 10 s , ae ~ 8 x 10 5 , and we 



(38) 



have a small parameter for our task 
T o(z) 

Here T c is the temperature of the solar center. For the 
whole Sun this parameter e is changed in the interval 
10 -8 < e < 10 -5 . e characterizes the degree of nonadia- 
baticity of the waves, it is defined as the ratio of the wave 
period to the reciprocal of the Kclvin-Helmholtz time. 

3. Adiabatic case 

For idealized adiabatic waves (e = 0) we have a second 
order equation, 

6" + fa®' + fa® = 0. (39) 
Introducing a new variable it may be written in standard 
form 



e = y^o /2 «rf 2 , y"{z) + i(z)Y = 0, 

where 



fa 



1 d 



(40) 



(41) 



fa and a are defined by Eqs. (37') and d = RqN 2 /g. The 
behavior of the function I(z) gives a possibility to ana- 
lyze qualitatively the waves in the solar interior and the 
boundaries from which the waves are reflected and become 
trapped. If / > we have oscillating solutions and if I < 
the waves are exponentially decreasing (evanescent) with 
z. Now we shall breafly discuss the incompressible and 
compressible cases. 

As we consider adiabatic waves, the transition to the 
incompressible case can be done by c 2 — > 00. In this case 



a3 p — ► 0, d — > — ae p , |A| — > 00 and instead of Eq. (py) we 



p 

obtain 



I 



5 1/5 

1 ) --ee T se Q -- I -ee T 



aen ,(42) 



where fc a = k x RQ and aen = dlnft/dz. Here there are three 
possibilities: 

1) absolute rigid body rotation, j3 — aen = 0. In this case 
/ < and we have no cavity solutions; 

2) rigid rotation with respect to latitude (f3 — 0), but 
'vertical' differential rotation with respect to the radius 
(aen 7^ 0). It follows from Eq. ( |42| ) that in this case for 
oscillating solutions we must have aen > as ae^ < 
is obeyed in the inner part of the Sun. It means that 
solar rotation must have a decreasing speed toward the 
center (dft/dz > 0). Then for / > the condition 
10|aeT|aen > (5aer/2 — aen) 2 must be fulfilled; 

3) The most realistic case for the Sun is a differential ro- 
tation in both directions ((3^0 and dft/dz ^ 0). In this 
case we have various chances to get trapped waves. We are 



interested in oscillations with 1 



10" 



x s 



(years) and 



k x w 10 2 (A » 3 x 10 4 km, e.g. sunspots). For solar condi- 
tions (ft w 2.86 x 10~ 6 s" 1 and A£ ax » 6 x 10~ 6 s" 2 ) we 
have k x u/2Q, « 1 and A 2 lax /w 2 w 10 10 . So the dominant 
term in Eq. (^) is the first one. To have an oscillating so- 
lution (I > 0) it is sufficient to have a very slow latitudinal 
differential rotation, (3 > 10 -8 w O(oj [s™ 1 ]). 

For waves running in opposite direction to rotation in 
the azimuth (f3k x > 0), the cavities (trapped wave area) 
can form between the bottom of the convective zone and 
almost the center of the Sun as well as in the outer part 
of the Sun where A 2 > 0. 

Waves propagating parallel to rotation may be trapped 
only in the convective zone (A 2 < 0). To the outer and 
inner sides from the convective zone the amplitude of these 
waves decrease exponentially. 

For the incompressible case it is easy to solve the 
eigenvalue problem of the cavity oscillations, because 
Eq. (42) has no singularity. Such a task has been solved 
by Oraevsky & Dzhalilov (1997). However, in the nonadi- 
abatic case we cannot apply the limit c 2 — > 00. Therefore 
we have to investigate the more complicated compressible 
case. 

To investigate the function I(z) given by Eq. (^l|) in a 
compressible plasma we need the orders of the quantitities 
entering the function I(z). To estimate these values let us 
consider a linear profile of temperature, To « T c (l — /3tz), 
where the gradient j3x = 1 — T c g/T c ps 1, T c ff is the 
effective temperature. Then we have a limit for the pa- 
rameter 86t from the center (z — 0) up to the sur- 
face (z = 1): 1 < — sct < 10 3 . The other parameters 
have the same order, ae PiP ~ 0(bst)- Then we get also: 



-ae 2 -' ae^ ~ ae P aeT, se' p ~ 3eT(ae p + asy). Now we 



can estimate the sign behavior of I(z). We shall consider 
more characteristic places of the Sun. In the following the 
condition /3k x > will be supposed. 
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At the center, where z — > 0, g — > 
0, c « 350km/s, if A > 1 we have 

7 « — pk x — sen —sex — aeo 

2ft 2 4 I 2 



0, iV 2 -> 0, 38 p 



1 2 

2^ 



< 0. 



In the middle part of the Sun (between the edge of 
the core and the bottom of the convective zone), where 
iV 2 w iV^ ax and ae p w const w — 10, the dominant term in 
the function I(z) is the first term of function ipQ. Hence 
we have I > such as in the incompressible case . 

The area around the convective zone is more compli- 
cated. The function I{z) has singularities at the points 
where d = and a = 0; d is connected with the Brunt- 
Vaisala frequency and a is defined mainly by the rotation 
gradient, Eqs. (37'). The two points of d = correspond 
to the bottom and the upper boundary of the convective 
zone. The function a{z) is defined mainly by two terms 
7 - 1 d' 



-A- 



(43) 



7 d 2 ' 

It is clear that as d! < at the lower boundary of the 
convective zone, the two zeros of a(z) are located around 
this boundary where d = 0. We can easily conclude that 

2 



lim I(z) 

a^0 



lim I(z) « 2 I — 
d^o V d 



a 
r\ 2 



(44) 



In the convective zone (A^ 2 < 0) we have I < 0. At 
the surface of the Sun we have again A^ 2 > uj 2 and hence 
I > 0. 

Thus if we approach the convective zone from below 
there exists a sequence of special points: four times the 
function I(z) crosses zero (turning points), and between 
these zeros the singular points are placed (see Fig.l). The 
turning points are determined practically by the condition 
of N 2 (z) — uj 2 = and the singularities by the conditions 
a = (circles in Fig.l) and N 2 = (boxes in Fig.l). Due 
to the very low frequency and the sharp decrease of A~ 2 (z) 
the turning point of the main inner cavity zt is very close 
to the first singular point where a = 0. 

In this way for waves with k x > a main large in- 
ternal cavity is placed practically between the center and 
the bottom of the tachocline. The solar atmosphere is a 
wave-propagating zone. Between the inner cavity and the 
solar surface a dark convective tunnel is placed; a very nar- 
row wave-trapping zone around the bottom of convective 
zone is also possible. It is clear that a tunneling of waves 
across the magnetized and turbulent convective area to the 
surface is probably possible. For waves with k x < the 
convective zone becomes a cavity. In this case the waves 
cannot be propagating at the solar surface. 

Not all singularities of the wave potential I(z) are sin- 
gular levels of the physical variables. Around the point Zd 
where d(zd) = we can write d fts d 1 '(z — Zd). Then we have 
the equation x 2 Y" + 2Y = 0, where x = z — Zd- The solu- 
tions of this equation are Y± t 2 — y/xx ±l ^^ 2 . As a ~ dr 2 
we get from Eqs. (|| + |c]) 'that p'/p ~ 6 ~ Y lj2 -> 



surface 




Fig. 1. Scheme of the dependence of the wave potential 
I(z) defined by Eq. (^|) on the distance normalized to 
the solar radius. Zeros of this function I(z) = (turn- 
ing points) divide the wave zone into transparent (cavity) 
and opaque (tunnel) parts. In the upper part of the Sun 
there are singular points of the function I(z) which are lo- 
cated between the turning points: the circles correspond to 
a(z) = and the boxes correspond to N 2 (z) = (bound- 
aries of the convective zone) . The narrow area at the bot- 
tom of the convective zone between circles and boxes is the 
tachocline. The main internal cavity comprises the whole 
radiative interior as the convective zone becomes opaque 
for the waves. 



if z — > Zd- It means, that the boundaries of the convec- 
tive zone are not singular levels for the initial physical 
variables. 

Another situation exists at the point z — z a where a = 
0. If we denote now x as x — z — z a , our Eq. (^o|) around 
the point x w is x 2 Y" ~ Y"3/4 = 0, the solutions of which 
are Y\ — x 3 ^ 2 , Y 2 = x~ x l 2 . Hence, for these solutions we 
have 0i x 2 7 2 ~ const. Then we get from Eq. (33) 
that the second solution in p'/po diverges at a = 0. 



There exist methods to construct asymptotic solutions 
of differential equations of second order with a singu- 
lar turning point. However, we are now returning to our 
fourth order Eq. (|3^) for two reasons: our intent is to 
consider the instability problem of the eigenmodes, and 
consequently in the complex ui plane the singularity at 
a(z,to) = is removed from the real z— axis. 

4. Asymptotic solutions 

The existence of the small parameter e in Eq. ( |37j ) allows 
us to apply asymptotic methods to solve this equation. 
Here we shall construct the inner cavity solutions only. 



As it has been discussed above the coefficients of Eq. (37) 
vary over a wide range. Very crudely, we have the fol- 
lowing estimates from the center of the Sun to the bot- 
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torn of the convective zone : e w O(10 8 — 10 5 ), sbt ~ 
O(l-10 3 ), tp 3 « 0(ae T ), <^ 2 » O(ae^), p x r> 0(ae|), » 
O(ae^), -01 ~ 0(a&r), and 0o ~ 0(ae§,). Assuming that 
eas^ w O(10™ 8 — 10) we can separate the variable part of 
e as e — £T c /Tq(z) and rewrite Eq. ( |37j ) in a convenient 
form: 

e(9"" + $ 3 e'") + $ 2 e" + *i0' + $ o = 0, (45) 
T T 

T 

$2 = =-(l+E^ 2 ), $3 = ¥3- 
1 c 

The general solution of this equation describes the slow 
and fast oscillating motions coupled with each other 
(quasi-adiabatic and dissipative modes). Every of these 
solutions is searched in a different way. 



Jf..l. Dissipative modes 

The solutions corresponding to the dissipative modes are 
searched for in the form: 

e = B (z) (l + £^ 2 B 1 +£B 2 + ---)e i " /2 r s ^ dz , (46) 

where S(z) and .80,1,2,... i z ) are unknown regular functions. 
Inserting the formal solution ( ^6|) into Eq. ([l5]) we may 
find all the unknown functions by the usual methods. In 
particular we have 



$2 



81 = - 



Fn = 



SB = 

Fx 
SB $ 2 

01 



Q-3/* e ifF d Z 



dz, 



$3 = 

$2 



■ £(f! 



<P3, 



(47) 
(48) 

(49) 



1 + £<fi2 

Fi = S(2SB'^ + 3S'B' + S"B ) + 3[S(SB )']' + 

+3$ 3 S(SB )' + $!8 + $ B - (50) 
In the derivation of these functions we have used the ob- 
vious condition £~ 1 / 2 SBq ^> B' . It may easily be shown 
that at the boundaries of the convective zone where d = 
the solution (Q) is limited. To solve the eigenvalue prob- 
lem we can restrict ourselves to the main branch of the 
asymptotic solutions of (|46|): 







1 



5/4 



i f F dz ±i f°(l +lp2 )i/ 2 dz 



(51) 



T H-e<^ 2/ 

In order to get the final form of this solution let us consider 
the extreme case e — > 0. Then we have 

exp \±i / -= , (52) 



where 1/e = i/\e\ and 1/v 7 ^ = ±(1 + i)/\/2]e|. It is clear 
that the turning points for these waves do not exist in the 
usual sence. It means, that these waves must propagate 
from a source down and up by decreasing their ampli- 
tudes. Hence only the 'transition' turning point for the 
dissipative modes can appear. The details of this question 
have been discussed by us in a recent paper (Dzhalilov et 
al. 2000) for the p-modes leakage problem. The place of the 
location of the sources of these waves can be taken in those 
places where the adiabatic approximation fails. Around 



the upper turning point z — z t the waves are rather nona- 
diabatic, where their vertical wavelengh is large. 

Thus we have a radiation boundary condition for the 
dissipative modes: from the point z — z t (turning point of 
adiabatic waves) the dissipative waves are radiated. In all 
directions the amplitudes of these modes must decrease, 
and while setting e = these modes must disappear. The 
following solution obeys these conditions: 



= 07 



T, 



1/4 



O j-Jo 



Jo(z) 



J(z) 



J(z) 



(l + £<p 2 )^ 

£G(z)dz, G(z) 
1/2 

¥>2 



Pi - 0112 
1 + £(/?2 



(53) 
(54) 



1 



dz, if z > Zt, 



1 

- + <P2 



1/2 



dz, if z < Zt- 



(55) 



is choosen, Re( J) < 



Here the branch with V2i = +(1 + i) 
0, and Tq(z) is normalized to T c . 

4-2. Quasi-adiabatic modes 

Now we will search for 'slow' quasi-adiabatic wave solu- 
tions in the form 







5> 

n=0 



Yn(z), 



(56) 



where the functions Y n (z) are supposed to be smooth. 
Inserting this formal solution into Eq. (45) we receive the 
reccurent differential relations between the Y n functions: 
*2^ " + + *o*o = 0, (57) 

d \ d 3 



dz 



-^-i, (58) 



where n = 1, 2, • • •. For our aim we are interested in the 
main value of the solution (J56l) only. Therefore for small 



£ Eq. (57) describes the main properties of the quasi- 



adiabatic waves. To obtain the standard form of this equa- 
tion we introduce a new variable 



— -dz 

$2 



W = T$ / *-/a\d\e Jo W, (59) 



Then we have 



Yo = cxp f-i 

where Jo(z) is defined by Eq. (fe 
W" + (\ 2 p + <p)W = 0, (60) 
where A 2 = c 2 A/c 2 , ujq = luRq/cq, cq is given for the 
temperature in the solar center, and 
A - a n u 2 ( N 2 



OtQ = 1 



(p{z) 



l + £f 2 A 

dz J d 
bi (ai + 35t + b[/bi) 



<v,„ 



1 + efc 2 



- q 
w ~d 



(61) 
(62) 



■ £ip 



1 + £<p 2 

-l^ 1+ eG) 2 -~(^+£G). (63) 

Here A =const is a large spectral parameter, and (p{z) 
is a smooth function. The asymptotic theory for differ- 
ential equations of second order such as Eq. ( |60| ) is well 
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developed for real p. In our case the function p is complex- 
valued with a complex spectral parameter u>. In such a sit- 
uation the direct application of the standard asymptotic 



theory to Eq. (60) is impossible. For the complex- valued p 
the turning points of Eq. ([30|) are shifted from those of the 
adiabatic Eq. ([l0|). Such small shifts are very important 
for the solution of the frequency discrepancy problem of 
helioseismology (Dzhalilov et al. 2000). For our task here 
the small shifts of the turning points are unimportant. 
Thus we take for Eq. (j^Tj) the same turning points which 
we have discussed in the section of adiabatic waves. For 
the solutions we apply the same method as earlier for the 
p-modes (Dzhalilov et al. 2000). We get the solutions 

1/4 

ir=[-| [dAi(-0 + CM-O] , z< zt, (64) 



t = 



2" 



2/3 



v = x 



where £,'( z t) ^ and 

P 



£ ^pdz, i\z) = -aJ|, (65) 



W 



[Ci^O+Mfc)], z>z t , 



(66) 



t = 



2" 



2/3 



v = x 




Here the regular functions and -Bj(£) are the Airy 

functions of first and second kind and C1.2 =const. In the 
derivation of the solutions (^34| + ^36|) the branch — 1 = 
exp(i7r) has been taken. 

Now we have to determine from the boundary con- 
ditions one of the unknown constants C\^. Here we as- 
sume the next simplification. In the real solar situation we 
should include the tunneling of waves through the opaque 
convective zone to the transparent solar surface. However, 
in this work we do not complicate the situation by includ- 
ing this important effect. Here we are interested in the 
eigenoscillation spectrum of the main interior cavity. So 
our solutions must be finite in the whole domain of inte- 
gration: < z < 00. In the limit of z — > 00 we have the 
asymptotics 

MO » J-r^e"", « 4=r 1/4 e", (67) 



2^~ 

where 77 is complex and Re(?y) — * 00. From here we have 
the condition C2 = 0, and the solutions (|5| + ^36|) are ex- 
pressed only by the function For small e we have 
the cavity solution for the quasi-adiabatic waves from 
Eqs. (H), ©, and 



6 = e NA ^T* /4 Va~\d\ 



1/4 



(68) 



where T$(z) is normalized to T c . This solution is limited 
at the boundaries of the convective zone (d = 0). Now 
we can define the general solution of our general Eq. ( |37| ) 
as a superposition of the solutions for the non-adiabatic 
(Eq. (|6"8|)) and the dissipative (Eq. (|53|)) cases: 

e = CxQ NA + c 2 e D , (69) 

where the new arbitrary constants C\ and C2 must be 
determined from the boundary conditions. 



5. Boundary value problem 

In this section we impose physically reasonable boundary 
conditions to the general solution of (jf^]) to determine the 
spectrum of the eigenoscillations of the inner cavity. 



5.1. Solar center 

At the solar center (z — 0) where g — > 0, ee p — > 0, 
and asx.p = 0(1) we apply a rigid boundary condition: 
v z (0) = 0- As the function Q(z) is finite at the center 
Eq. d|) reads as d(p'/p )/dz = 0. Using Eqs. (§| + £l|) 
with dimensionless parameters wc have 

0. (70) 



dz " 



e 



It is easy to show that for both solutions the conditions 
Q% a /&na ~const, &'[)/& d ~ const are obeyed. Then 
the condition ([70|) is changed to 

dz z =o 

Now we insert the general solution ( |69| ) into this condition 
and receive the ratio between the unknown coefficients: 



9i 



®'na(0) 



(72) 



To obtain the dispersion relation we need one more bound- 
ary condition. 

5.2. Bottom of the convective zone 

Let us consider the location of the upper turning point 



(z = Zt): there the equation N 2 



is fulfilled, from 



which we may define the parameter d, d rj 0(RqLu 2 / g) w 
10 -9 w 0. Here we can apply a free boundary condition, 
because at the bottom of the convective zone the Brunt- 
Vaisala, frequency and the rotation rate fl(z, 9) change 
very sharp. (A more realistic approach would be to in- 
clude the tunneling of the waves.) At the free surface the 
Lagrangian pressure must be constant, dp/dt — 0, where 
p = po + p' is the total pressure. From here we have 



iwc 2 p' 



z t - 

9 Po 
From Eqs. (|| 
in dimensionless parameters 
qui 2 ( N 2 



(73) 



27) we define the next formulae written 



^ = e 

Po 

1 

+ d 



• e 



— 33 



dz 2 
d 
dz 



e 



dz 2 



Ik 
c 



fid, 



1 



dz 2 



e 



p 

nq — 

Po 



(75) 



Having in mind that at z Zt we have a —qd'/d 2 and 
d! < 0; neglecting small values we get from Eqs. (f73|) + 
75|) 

d 2 



1 



dz 2 



O = nq 



P_ 

'Po' 



(76) 
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Now excluding p' /po from here and from Eq. ( J74| ) we get 
the second boundary condition 

z—z t 

Inserting the general solution (p9) into this condition we 
get 

C2 6' NA (zt) (78) 



5.3. Dispersion relation 



we get the general dispersion relation 



From Eqs. (|7j 

e^(o) = ei 3 (o)^M^. (79) 

To simplify this equation the asymptotic expansions of the 
Airy functions can be used. At the center of the Sun the 
arguments of the Airy functions become large: |£| >> 1. 
Then the main values of the asymptotic expansions yield 
(Eqs. 68 + 53) 

i\ 2 p, 



■ cos 



+ 1) 



e' D (o) ~ -—exp^y^ ^d, 

where 770 = ??(0) (Eq. 65) and p c = p(0) 
At the turning point z — z t we have 

3/4 



(81) 

-^0 (Eq- si). 

(-A) 1/2 (p'(^)) 1/6 A 5/6 4(0), (82) 

(83) 



T 

,T C / ye 

where T t = T {z t ). Now setting Eqs. (^0[) -(|§) into @ 
we get 

TT, . , ,p'(z t )M'9) 1/2 



cos(r?o + -) 



A 7 /6p, 



1/4 



x 



x exp 




2|£| JO 

The right-hand side of this equation is negligible small as 
Ai'(0) = -0.259, A > 1 (Eq. 60), q « (7- l)/ 7 (Eq. 37'), 
d'(z t ) w 16, and mainly e » 10~ 8 (Eq. 38). Thus we have 
finally using Eq. (65) 

xj^ Vpjz)dz = 7r(n-^). (85) 

We may rewrite this dispersion relation for the i?-modes 
in a visually more convenient form 

Ph (86) 



I(oj,k x 



A{z,oj) = 



N 



1 - a /A 



I A 



7p 



dz, 



7r(n-3/4) 



1 



^2 



Here n — 1,2,- •• are the radial node numbers, for the 
other parameters see Eqs. (37' + 62). In the adiabatic limit 
(e = and then Im(tj) = 0) for very small frequencies 



(ui 2 <C N 2 ) and for incompressible motions (c 2 — > 00) 
we get the classical dispersion relation for Rossby waves: 
W ~ fik x j{k\ + k 2 ) with fc 2 = cr 2 . 

6. Instability of 22-year, 4000-yr, and 
quasi-biennial oscillations 

Now we can determine the complex eigenfrequencies 
lu = u> n (/3, k x , ky, a) from the integral dispersion relation 
Eq. @ for the given fl Q = 2.86 x lO^s" 1 . Remind 
that (3 is the latitudinal gradient of the rotation rate, 
f2 ps Sl(l — Ph/Rq), and a determines the power exponent 
of the temperature of the nuclear reactions, Q ~ p 2 T a 
(Eq. 5). For the p-p reactions we will use a — 4. The 
dependence of this equation on the free parameter k y is 
a drawback of this equation. In real situations the waves 
cannot be progressive across the shear in j/-direction. In 
the general case the complex k y = k y (k x , (3, ui) is a solution 
of the two-dimensional eigenvalue problem. To simplify 
our task for a very small latitudinal gradient of rotation 



(80) (3 we have introduced the free parameter k y . To simulate 
the decay of the wave amplitudes towards the pole, in y 
direction, we shall consider only the case with real k 2 < . 
A negative spectral parameter k 2 could also be confirmed 
by the tidal equation of Laplace (e.g. Lee & Saio 1997) for 
very low frequencies. This equation includes the influence 
of sphericity on the angular dependence of the eigenfunc- 
tions, if simple rigid rotation is considered in the frame 
of the traditional geophysical approximation (which fails 
near the equator). 

Ando (1985) has derived a local dispersion relation for 
waves around the stellar equator. Solving his Eq. (15) with 
respect to k (that is our k y ) for low frequencies (uj/20, <§; 



1) we get two solutions: k 2 = —m 2 and k 2 — ~ m 2 N 2 /4£l 2 
where m is the azimuthal wave number. As iV 2 /4J7 2 ^> 1, 
the second solution is strongly damped in y— direction. 



(84) Our case corresponds to the first solution, k = 



From geophysical applications (Pedlosky 1982; Gill 
1982) we infer that k 2 (k x ) ~ —k 2 is obeyed for the 
ageostrophic wave propagation across the shear flow. Re- 
ally, it is seen from Fig. 2, that the mode separation with 
respect to n is essential for the solar situation if the con- 
dition k 2 w — k 2 is fulfilled. In Fig. 2 the Re(uj(k y )) depen- 
dence is shown for the case k x = 1 and (3 = 10~ 6 for ex- 
ample. Here, a very small imaginary part of the frequency 
is supposed, fj< 1, where 

uj = uj r (l + in), 77 = Im(w)/Re(w). (87) 
In the case k 2 /k 2 7^ —1 we have a continuous spectrum 
which is physically not interesting. So during this work 
we have excluded the parameter k 2 from our dispersion 
relation (^ using k 2 = — k 2 . As a result the frequencies 
depend only on the product (3k x ; the dependence on a is 
very weak. 

It follows from Eq. ( |86| ) that for a given (3k x the fre- 
quencies u) n decrease with increasing harmonic number n 
(u> ~ 1/n 2 ) in the adiabatic case (e = 0). At very high 
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Fig. 2. Mode separation region versus k y = fc a i? Q . On the 
vertical axis the absolute values of the frequencies normal- 
ized to LO22 are shown. The maxima from top to down in 
the curves correspond to n=l, 3, 5, 19. Calculations 
are done for the case k x = 1, (3 = 1CP 6 , rj = 1CP 5 . This 
figure does not yet show the eigenfrequencies, but the do- 
main of ky where the n-dependence of the frequencies is 
obvious. This domain is k 2 — k 2 ,. 



tl — > oo we have almost steady motions with a frequency 
uj — > 0. Howewer, the situation is changed if dissipation 
is taken into account, e ^ 0. In this case small-scale mo- 
tions are quickly damped, very low frequency oscillations 
with high n could not become trapped and cannot mani- 
fest themselves as eigenmodes. Really, as the nonadiabatic 
parameter e ~ at very low frequencies the solutions 
of the initial equations must have a dissipative character. 
An investigation of the initial equations for a steady flow 
with d/dt = is a separate task. Here we restrict ourselves 
only to our dispersion relation Eq. ( |8q ) with increasing n. 
We can rewrite this equation as lfnr » lu/I 2 . With de- 
creasing uj the right-hand part uj// 2 decreases, and then it 
becomes independent from uj. It means, that n is limited: 
n < rimax- Hence, it is a results of the dissipative effects, 
that only limited modes become trapped. The value of 
rimax depends on k x and j3. In Fig. 3 the n miac (k x ,P) de- 
pendence is presented. It follows from this figure that for a 
given k x the number n max increases strongly by increasing 
the rotation gradient (3. So, if the solar interior is rotating 
similar to a solid-body, very long-period oscillations (al- 
most steady flows) should be suppressed. As the values of 
n max are sufficiently high, the accuracy of the asymptotic 
solutions is high. 

Now we consider which spectrum of trapped waves 
with 1 < n < ri max is possible. To calculate the wave 
spectrum from Eq. (§3) we use the standard model of 
Stix (private communication; Stix & Slaley 1990) for the 
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Fig. 3. The harmonic numbers (n max ) limited by the 
nonadiabatic effects versus k x and (3. The quantities on 
the curves correspond to the gradients of the rotation rate 
(3. Such a limit for high n does not exist for the adiabatic 
waves. 

solar interior. For the special case k x = k x R Q — 100 (i.e. 
X x 1=3 40 000 km) the n-dependence of the real and imagi- 
nary parts of the eigenfrequencies are shown in Figs. 4a, b. 
The calculations were done for different small values of the 
rotation rate gradient /3, covering a wide range: j3 = 10~ 7 - 
10 -4 . This is done because we know from helioseismology 
only that the (3 parameter is small, but the exact value is 
not yet known. In Figs. 4a, b uj r =Re(w) and uji =lm(uj) are 
normalized to the cycle frequency of the 22-year oscilla- 
tions: uj-2.2 = 0.91 x 10 _8 s _1 = 27r(1.45nHz). As expected 
the frequencies decrease with n and increase with (3 or k x . 
The imaginary parts Wi oscillate around the zero value: 
if u>i > 0, the waves are unstable and their amplitudes 
increase with time; in the opposite case, if a/j < 0, we 
have stable/damped waves. In Fig. 4b we have two posi- 
tive maxima: the first corresponds to short-period oscilla- 
tions of 1-3 yr ('quasi-biennial modes') and the second one 
to medium-period oscillations of 18-30-yr ('22-yr modes'). 
The position of the quasi-biennial modes versus n is stable 
and is ~ 15. For the 22-yr modes n.22 is slightly increased 
with an increase of (3. It is seen from Fig. 4a that n must 
increase to keep the same frequency with increasing (3. For 
smaller frequencies this shift is stronger. Fig. 4b shows that 
the instability gets stronger if (3 increases: unstable waves 
become more unstable and damped waves are stronger 
damped. Waves with high n > 200 at very low frequencies 
also show instability which cannot be shown in Fig. 4b due 
to the scaling. 

Of course, f3 will change with the radius in the real 
solar radiative interior. Hence, those places, where (3 be- 
comes relatively large, may become sources of unstable 
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Fig. 4. Real (uj r ) and imaginary (uii) parts of the eigen- 
frequencies ((a) and (b), respectively) normalized to u>22 
versus the radial harmonic number n. The numbers on 
the curves 1, 2, ... , 10 correspond to different values 
of the latitudinal gradient of the rotation rate /?i-io = 
10" 4 ,7 x 1(T 5 ,5 x 1CT 5 ,3 x lCT 5 ,l(r 5 ,7 x 1(T 6 , 5 x 
1(T 6 , 3 x 1CT 6 , 1(T 6 , 1CT 7 . Positive uj 1 corresponds to the 
growth rate of the waves due to the £-mechanism. For neg- 
ative uji the waves become damped in time due to radiative 
losses. 

waves. Our calculations can easily be generalized for any 
k x as Eq. ( p6| ) is a function of (3k x only. 

To characterize the mode instability, the behavior of 
the parameter r\ is more important, rj is the growth rate 
(increment) of the instability of the modes if rj > and 
the damping rate (decrement) if rj < 0. In Fig. 5 we present 
the rj(ui r ) dependence for the whole range of frequencies 
(uj r < Qq) for which our asymptotic theory is valid. We 
have three distinguished global maxima for the growth 




100.0000 



Fig. 5. Growth rate of the instabilities of the eigenmodes 
of the differentially rotating Sun. Three modes with pe- 
riods of ss 1-3 yr, 18-30 yr (these ranges depend on (3k x ) 
with a small additional peak at 100 yr, and 1500-20 000 yr 
(independent of pk x ) become maximum unstable (rj > 0) 
for high orders n. The solid, dashed, and dotted lines cor- 
respond to the small latitudinal gradients of the rotation 
rate of (3 = 7 x 10~ 6 , 5 x 10" 6 , and 10~ 5 at k x = 100. 



rate rj > which correspond to period ranges of « 1- 
3yr, 18-30 yr with a small additional peak at 100 yr, and 
1500-20 000 yr ('4000-yr modes') of the eigenmodes. These 
modes have radial node numbers near n = 850, 60, and 
15, respectively. It was already mentioned, that for the 
unstable modes the value of n is slightly changed with a 
change of (3 . So for the unstable modes n is high and 
the asymptotic results are reliable. The growth rates of 
the 22-yr oscillations are always greater than those of the 
quasi-biennial modes. The characteristic growing time for 
the 22-year modes is « 1000 yr as r\ ~ 0.003. 

Now we consider in which range of n and j3k x the 
unstable modes are located. For this aim we fix the fre- 
quency as uo r = LO22 = 0.91 x 10 _8 s _1 , and for a given 
n < n max we find the complex root of the dispersion re- 
lation Eq. d86| ) . The parameters j3k x and rj correspond to 
this root. The same calculations were done for the 2-yr 
oscillations. The results are shown in Fig. 6. Here the solid 
lines correspond to the 22-yr and the dashed lines to the 
quasi-biennial modes. A sharp increase of n to n max with 
/3k x in Fig. 6a indicates some upper limit for the param- 
eter (3k x . For the 22-yr and 2-yr modes these limits are 
approximately 3 x 10~ 3 and 3 x 10~ 2 , respectively. Fig. 6b 
(where the increment /decrement is presented) shows for 
which values of j3k x and n the 22-yr and 2-yr modes be- 
come unstable. Two maxima for each mode appear in rj. 
These are for the 22-year modes: 
n = 15, /3k x = 9 x 10" 5 , rj = 0.004 and 
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Fig. 6. The existence area of the 22-yr (solid lines) and 
2-yr (dashed lines) eigenmodes versus the parameters 
n, (3k Xl and r\. Fig. 6a shows that the wavenumber and the 
rotation gradients {(3k x ) are limited as n < n max . The two 
narrow maxima of the growth rate r\ in Fig. 6b indicate 
that for the strongly fixed values of f3k x the 22-yr and 
2-yr modes become unstable (see text). 



n = 55, (3k x = 7x 1(T 4 , r, = 0.003; 
for the 2-yr modes: 

n = 14, f3k x = 5x 10~ 4 , t] = 0.002 and 
n = 65, (3k x = 10~ 2 , i] = 0.001 . 

These estimates could be used to define a possible 
value of (3 inside of the Sun. This will be possible if we 
can identify these modes from observations. For instance, 
at k x 100 (sunspot scale) a gradient of the rotation rate 
with = 9 x 10~ 7 or {3 — 7 x 10~ 6 is needed to excite 
the 22-yr modes. Here the possibility of an excitation of 
both the n — 15 and the n = 55 modes is not excluded 



if (3 is changed with the radius. A better way to define j3 
would be to identify both the 22-yr and the quasi-biennial 
modes, with different k x at the surface of the Sun. 



7. Conclusions 

In the present paper we have shown that toroidal eddy 
flows which are degenerated in a non-rotating fluid can be- 
come a reservoir of various branches of oscillatory modes 
when the degeneracy is removed by rotation. The mech- 
anism depends on the condition for the existence and al- 
teration of the relative vorticity as well as on the stel- 
lar rotation rate and its gradient. Apparently at least 
for slowly rotating stars (£1 < il g ) the rotation waves 
could be divided into two types: r-modes with high fre- 
quencies (lu < Q) which are independent from the inner 
structure and mainly caused by geometrical effects, and 
the i?-modes with low frequencies (w <C il) which depend 
on the inner structure and are considered in the present 
paper. This classification is similar to that of /- and g- 
modes or to that of surface and body tube modes of mag- 
netic cylinders. Note that the properties of Coriolis forces 
and ponderomotive forces in MHD are very similar to each 
other. Both rotation modes are prototypes of the geophys- 
ical Rossby waves. 

We investigated the instability problem of the i?-modes 
sustained by a very small latitudinal gradient of the ro- 
tation rate in the solar radiative interior. Among the 
eigenoscillations three modes with periods of « 1-3 yr, 18- 
30 yr, and 1500-20 000 yr turn out to be maximum un- 
stable to the e-mechanism. Here the smoothing effect is 
the radiative damping. All of these instabilities are in the 
range of high radial node numbers n which indicates that 
the applicability of the asymptotic solution is satisfied. 

The 22-yr modes with a growing time of ~ 1000 yr are 
of particular interest with respect to the solar activity cy- 
cle problem. In the simpler case when adiabatic i?-modes 
are considered in an incompressible fluid, a n in Eq. (86) 
is independent of the wave number and of the frequency 
for very low frequencies. Then in the azimuthal direction 
the phase and group velocities are 

(3 ^_kl-kl-ol 

Vp X ~ ~2 ( j) i ^gx ~ P 



kl+ar " {k'i+alf 

respectively, where the velocities are normalized to VLqRq 
w 1.6km/s. In our case (3 > 0, i.e. the angular velocity 
is decreasing towards the pole, similar to the behavior at 
the solar surface. The x-axis in our coordinate system is 
directed opposite to the direction of rotation, v px > and 



< 0, moreover, kf. 



-k"i and for finite n we have k'i > 



<7„ (see Fig. 2). Then the retrograde R- modes transport 
energy along the direction of rotation. Our treatment in a 
Cartesian coordinate system does not allow to determine 
the direction of energy transport by wave packets relative 
to the equator in the meridian plane. The estimate of fc 2 , ~ 
— k x is crude, and to determine the exact dependence on 
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k y — k y (k x ,uj 7 fi) the 2D boundary value problem must be 
solved. 

A nice property of the Rossby waves is that every 
monochromatic mode is a solution of the full nonlinear 
hydrodynamic equations. It means, that we should expect 
the development of nonlinear i?-modes with large ampli- 
tudes. We could also expect that just in this nonlinear 
regime the toroidal magnetic flux will be lifted from the 
upper boundary of the cavity (the tachocline) to the sur- 
face. The energy release of the nonlinear waves could be 
accomplished by magnetic reconnection. Here it is pos- 
sible that toroidal currents are generated via a twist of 
toroidal magnetic field lines by the cyclonic flows of reg- 
ular i?-modes with fixed characteristics. Parker (1955) as 
well as Steenbeck et al. (1966; see also Krause & Radlcr 
1980) have suggested for the dynamo process that such a 
mechanism, the a— effect, is working by turbulent motions 
under the influence of Coriolis forces. 

Our present model points out the possibility of forced 
oscillations instead of a self-excited dynamo to solve the 
solar cycle problem, and this with the correct period of 22 
yr. Similar ideas are due to Tikhomolov (2001) who has 
recently suggested a hydrodynamic driving of the 11-yr 
sunspot cycle. We expect that in our model — contrary to 
classical dynamo models — a huge toroidal magnetic field 
of ss 10 5 G will no longer be required to explain the buoy- 
ant rise of magnetic flux tubes appearing at the surface 
with small tilt angles and at low latitudes: the external 
nonzero upflow produced by the regular vortical i?-modes 
could trigger the eruption of stable magnetic flux tubes 
stored in the overshoot region. There is still a smaller peak 
of the growth rate (Fig. 5) at 100 years; such a period is 
observed as a modulation of the 11/22 yr cycles. 

There is observational evidence for the short-period os- 
cillations as well: From helioseismic sounding Howe et al. 
(2000, 2001) have recently discovered variations of solar 
rotation with a period of 1.3 yr in the lower convective 
zone. Quasi- two year modes are very likely seen regularly 
in various solar data (e.g. Waldmeier 1973; Akioka et al. 
1987: Rivin & Obridko 1992). The existence of two mag- 
netic cycles (the main 22-yr and the quasi-biennial pe- 
riod) on the Sun has been reported by Benevolenskaya 
(1996; 1998). So far the origin of these modes was not 
yet clear. Terrestrial quasi-biennial oscillations have been 
clearly seen in tropical meteorological radiosonde data, 
and a possible solar origin by related phenomena in the 
solar interior, Rossby waves in particular, has been dis- 
cussed as well (Mclntire 1994). 

The long-period oscillations in the broad range 1.5 x 
10 3 -2.0 x 10 4 yr, with a maximum growth rate around 
4500 yr, could be the cause of abrupt changes of the global 
terrestrial climate in the past: Dansgaard-Oeschger events, 
these are abrupt onsets of warm periods during the last 
ice age, had mean distances of 4500 yr, but they were 
distributed over a larger period range, similar to that in 
our model, with shortest distances often around 1500 years 



(see, e.g., Ganopolski & Rahmstorf 2001). These events 
were caused by changes of the thermohaline circulation of 
the ocean, which in its turn were probably triggered by 
changes in the solar energy output. 
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